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The complex structures of warm and hot dense matter are essential to understand the behaviors of 
materials in high energy density physics processes and provide new features of matter constitutions. 
Here, along a new unified first-principle determined Hugoniot curve of iron from normal condensed 
condition up to 1 Gbar, the novel structures characterized by the ionic clusters and separated 
"electron bubbles" are revolutionarily unraveled using newly developed quantum Langevin molecular 
dynamics (QLMD). Subsistence of complex clusters, with bonds formed by inner shell electrons of 
neighbor ions, can persist in the time length of 50 femto-seconds dynamically with quantum flowing 
bubbles, which are produced by the interplay of Fermi electron degeneracy, the ionic coupling and 
the dynamical nature. With the inclusion of those complicated features in QLMD, the present data 
could serve as a first-principle benchmark in a wide range of temperatures and densities. 

PACS numbers: 52.65.Yy 52.27.Gr 62.50.-p 64.30.Ef 



The thermodynamic and structural properties of mat- 
ters at extreme conditions, so-called warm dense matter 
(WDM) and hot dense matter (EDM) in the field of high 
energy density physics (HEDP) [1], are both experimen- 
tal and theoretical challenges and critical to the com- 
prehension of the evolution and the internal structures 
of giant planets, stars, inertial confinement fusion (ICF) 
target capsule, and material science [1-7]. New physics 
discovered by high-power laser facilities such as National 
Ignition Facility (NIF) and nuclear fusion [8-13] requires 
understanding beyond the traditional condensed matter 
and atom (or plasma) physics [1, 6]. Recent laser-driven 
dynamical experiments and related theories show the 
existence of ordered electron-ion structures [13-20] and 
electronic bonds [11, 21] in WDM and HDM with x-ray 
methods, suggesting the need to take into account of the 
dynamics of local chemical environments. Meanwhile, 
studies of static high-pressure theories and experiments 
have led to the finding of " electron blobs" formed by va- 
lence electrons and new-type electronic bonds assisted 
by inner-shell electrons in cold aluminum and sodium at 
high pressures [22, 23]. However, for the lack of effective 
methods, few theoretical studies are carried out on the 
structures of complex materials in HEDP field, which are 
crucial for determining their physical properties such as 
energies, pressures and transport behaviors. 

In the HDM, the densities are comparable to or even 
much higher than the states of static compression and 
WDM, and temperatures are up to hundreds of eV [1]. 
The higher density effect could induce new features that 
cannot be seen in the normal WDM, while the higher 
temperature causes dynamical changes of the ionic con- 
figurations accompanied with dynamical electronic dis- 
tributions, which do not exist in the static highly com- 
pressed cold matters. The properties of these kinds of 



matters, including the equation of states (EOS), elec- 
tronic and ionic conductivities, viscosities and diffusions, 
and optical properties, are definitely dependent on the 
unknown details of the electron-ion structures under 
these extreme conditions. 

In order to shed light on the hidden features and 
controversial intrinsic dynamics from WDM to HDM, 
the electron-ion structures are calculated along the orig- 
inally determined principal Hugoniot curve of Fe us- 
ing newly developed first principles method-quantum 
Langevin molecular dynamics (QLMD) [24, 25]. For Fe, 
as one of the most abundant elements in universe and 
a typical complicated transitional metal, it is a long- 
standing challenge [3-5, 13, 26-30] to obtain the phys- 
ical properties such as EOS and electron-ion structures, 
because of the strong ionic coupHng, and high electronic 
degeneracy in a wide range of temperature and density. 
To date, previous experiments and statistical framework 
based theories generate abundant results with large di- 
vergence and uncertainty in EOS [26-30]. First princi- 
ples studies on the EOS and electronic properties of crys- 
talline or liquid Fe at high pressure and zero or relatively 
low temperature have been reported widely [3-5, 31, 32], 
and a few fixed density-temperature points on the Hugo- 
niot curve picked up from SESAME table were calcu- 
lated by quantum molecular dynamics (QMD) [25, 26]. 
However, none Hugoniot data beyond WDM from first 
principles are covered. 

Hugoniot curves are determined by both the pressures 
and the internal energies. QLMD or QMD, based on 
finite-temperature density functional theory (DFT) [33], 
can naturally include the effects of degeneracy and cou- 
pling contributing to the pressure and energy, and has 
been successfully applied to derive the EOS and dynami- 
cal properties of dense matter including Fe [2, 31, 32, 34- 



37]. Advantageously, QLMD, adopted in the Quantum 
Espresso package [38], can be extended to the HEDP 
field within the framework of ab initio by introducing 
electron- ion collision induced friction (EI-CIF) [25]. It 
is thus possible to accurately explore the details of the 
electron-ion structures in WDM and HDM as a power- 
ful tool. In the present work, 54 atoms are included in 
the super cell with 3x3x3 k-points below 10 eV and F 
point only at higher temperatures for the representation 
of Brillouin zone. Pseudopential (PP) with 16 electrons 
in the valence within the generalized-gradient approxi- 
mation (GGA) [39] is used, since the ionization degrees 
are less than 16 below 100 eV [25]. During the MD pro- 
cesses, the time steps are from 1 fs to 0.25 fs with in- 
creasing the temperature, and 2 ps time length is chosen 
to achieve the thermal stability state. After the thermal- 
ization, more than 2 ps time length is used to acquire 
the thermal properties such as EOS and ionic structures. 
Over 300 density-temperature points are calculated in 
order to get the EOS data. 

The EOSs of the temperatures from 0.1 eV to 100 eV 
and pressures up to 1 Gbar (1 Gbar = 100 TPa) on 
the both sides of the Hugoniot curve are obtained. We 
firstly discuss the electron-ion structures and the dynam- 
ics along this curve, and finally give the EOS data from 
normal condensed state to WDM and to HDM. One of 
the most spectacular physics here is the electronic struc- 
tures in WDM and HDM, which we knew few today. To 
dig out the undiscovered features, the electronic distribu- 
tions at the density-temperature point of (100 eV, 33.385 
g/cm^) on the calculated Hugoniot curve are displayed. 
It is shown in Fig. 1(a) that the formation of "blobs" of 
the valence electrons of cold iron at high pressure, which 
was proved in high density cold aluminum [22] . The con- 
traction at high pressure makes the ions close enough, 
with the free-like valence electrons occupying their left 
space. When the valence electrons become free, the in- 
ner 5, p-electrons will assist the bonding on Fe-Fe, as 
shown in the band structures in very recent high-density 
results [31]. This can be verified in the 2-dimensional 
density distribution in Fig. 1(b), where the valence elec- 
trons (low densities) are distributed between Fe ions, 
and the inner electrons form covalent bonds. How can 
this feature change when the dynamical effects are in- 
troduced? As shown in Fig. 1(c), the valence electron 
" blobs" tend to assemble together and form bigger " bub- 
bles" in the interspaces of Fe ions due to the temperature 
induced dynamics. These free electrons are not homoge- 
nous and behave as "quantum electron Hquids" flowing 
with ionic moving, which can not be described by the 
current atom and plasma models. It can be verified in 
the 2-dimensional picture in Fig. 1(d), where the free 
electrons (red color) are distributed in the interspaces 
of ions. Interestingly, there are clear-cut density over- 
laps among more than three ions at high density induced 
by the inner orbital electrons, indicating the existence of 
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FIG. 1: (Color) The electronic charge density (electron/A^) 
distributions of Iron, (a) and (b): the three and two dimen- 
sional in the (100) direction contour plot for the charge den- 
sity of Iron in bcc phase at (0 eV, 33.385 g/cm^); (c) and (d): 
the same contour plot of Iron at (100 eV, 33.385 g/cm^). 



many-body bonding formed by inner shell electrons. Fur- 
thermore, the covalent-like bonds are shown, which are 
likely formed by the 5 to p or sp to d electron transfer 
and their hybridization at high density [31]. 

The "electron bubbles" due to the interplay of 
the Fermi electron degeneracy, the ionic coupling and 
temperature-induced dynamics are obviously separated 
from the ions, which can be considered as two phases 
coexisting and flowing in HDM. This discovery is quite 
different from the physics in WDM and ideal plasma 
gases. In WDM, the electrons and ions are similar to the 
physics in condensed matter, where the electrons "move" 
clingingly to ions " simultaneously" and distribute around 
some specific atoms. In ideal plasma gases, the electrons 
and ions are all considered as classical non-interacting or 
weakly coupled particles. Here in dense and hot mat- 
ter, the ions can form networks distorting slowly in the 
quantum electron fluid, giving an alternative quantum 
mechanical case for the two phases. This feature with 
bonds should not be found in the simple elements such 
as hydrogen, whose electronic structures are much sim- 
pler, and there are no "inner-shell electrons" assisting 
the formation of new bonds and excluding the ionized 
valence electrons from the region of gathered nuclei due 
to the Fermi degeneracy. 

The electronic distributions show the complexity of the 
ionic structures, whose details and dynamics are still se- 
crets. Most importantly, the electronic structures are 
tightly dependent on the dynamics of ions and their col- 
lective behaviors. In order to understand the physics 
of structures, we select five temperature-density (T-D) 
points as shown in Fig. 2 along the new principal Hugo- 
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FIG. 2: (Color online) (a) RDF and (b) the distribution of ori- 
entation order parameters Q for the selected five temperature- 
density points along Hugoniot curve. The corresponding tem- 
peratures and densities are respectively (0.1, 1, 10, 40, and 
100) eV, and (10.1, 13.23, 18.75, 28.875, and 33.385) g/cm^ 

niot curve to uncover their hidden dynamic nature of 
the ionic structure. Their radial distribution functions 
(RDF) shown in Fig. 2(a) give evidences of the transition 
of ionic structures from long-range order to short-range 
one statistically. It is worth noting that even at T = 
100 eV, there is also one peak in RDF, indicating the 
existence of hidden ordered structures. 

Considering the short-range ordered structures at high 
temperatures, we borrow the language of the liquid struc- 
tures such as water and clusters to reveal the struc- 
tures in HEDP. The orientation order parameter Q = 
1 - i Z)j=i+i(cos% + is defined, where % is 

the angle formed by the lines of an ion and its nearest 
neighbors i and j (< 4). The value of Q varies from 
(in an ideal gas) to 1 (in a perfect tetrahedral net- 
work), which can be used as a measure of tetrahedrality 
for the local coordination structure [40, 41]. As shown in 
Fig. 2(b), the peak of the distribution of the parameter 
Q shifts from 0.45 to 0.35 with increasing temperature, 
indicating that the ionic structures with the ideal tetra- 
hedral network collapse indeed but there are still some 
similar topological structures at high temperatures. 

As for the topological structures, some ordered struc- 
tures survive even at 100 eV from the hints of Fig. 1 
and Fig. 2. However, there is not a clear minimum af- 
ter the first peak in the RDF in Fig. 2(a) except for 
0.1 eV. Therefore, a dissociation criterion based solely 
on a hard cutoff on Fe-Fe bond lengths would be op- 
tional. The probability distribution of coordination num- 
bers (CNs) can be good for analyzing the local geometry, 
which has been successfully used for estimating the dis- 
sociation for water molecules at high temperature [42]. 
Due to the short ordered structures in HDM [43], here 
we adopted the effective CNs (ECNs) idea[44, 45]. For 
low-symmetry structures where a particular atom is sur- 
rounded by atoms at different distances, the ECNs con- 
cept can be independent on the choice of the bond cutoff, 
and therefore provides a more accurate method to deter- 
mine possible structural trends in disordered structures. 
Thus, ECNs can lead to more clearly structural figures. 

For the above five T-D points, ECNs decrease from 
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FIG. 3: (Color online) The trajectories of ECNs for a specific 
atom during the time length of 0.5 ps simulations at different 
states of (a), 0.1 eV, (b), 1 eV, (c), 10 eV, and (d), 100 eV, 
respectively. 

10.99 to 2.93 A, whose corresponding average bond 
lengths are from 2.33 to 1.15 A. With increasing temper- 
atures and densities, ECNs decrease gradually, but retain 
larger than the unit even up to 100 eV. This fact indi- 
cates the existence of cluster-like or network-like struc- 
tures, which is consistent with the hints in Fig. 1. A 
part of ions catch only one nearest neighbor (ECN = 1) 
when temperature is high enough, indicating the forma- 
tion of two- ion chains. The distributed percentages of 
these "chains" are 13%, 26%, and 36% respectively for 
the temperatures of 10, 40 and 100 eV. Furthermore, the 
differences between the dynamical average bond lengths 
and ionic radii (from 2.10 to 1.41 A) show that the in- 
teratomic distances can not be simply described by the 
hard sphere model, and the average sphere space can not 
be only responsible for the density. 

The dynamics of these bond-network breaking and 
forming patterns, i.e., the dynamical dissociation of the 
Fe-Fe bonds, is also pivotal since an instantaneous topo- 
logical structure can not affect the observed physical 
properties obviously. In order to investigate the struc- 
tural dynamics, we firstly trace the topological networks 
around one specific atom, revealing the dynamic na- 
ture of the topological structures during the simulations, 
shown in Fig. 3. With the increasing temperatures, the 
movements of ionic positions in liquid states introduce 
some physical processes such as dissociations [42], re- 
sulting in the change of structures at different times. 
However, even though the structures exhibit dynami- 
cal changes, a considerable fraction of ions in warm and 
hot dense Fe with compact clusters (networks) can per- 
sist for a long time of a few tens of femto-seconds {fs) 
even at T = 100 eV. Because most values of ECN^ 
(ECNi represents the ECN of the ith atom) are located 
around ECN within the departure of 2, only three pat- 
terns of topological structures are statistically averaged: 
using ECN as the reference, for the ith atom, we assume 
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FIG. 4: (Color online) Principal Hugoniot of Fe, comparing 
the results from QLMD with other models and experiments. 
The corresponding temperatures are T = (0.1, 0.5, 1.0, 2.0, 
3.0, 5.0, 8.0, 10.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 50.0, 60.0, 
70.0, 80.0, 90.0, and 100.0) eV. 



ECNi=ECN-l when ECNi<ECN-l, and ECNi=ECN+l 
when ECNi> ECN+1, and ECNi=ECN when ECN- 
l<ECNi <ECN+1. Summing up the persisted time 
length for every structure, we can find that more than 
15% of the topological structures forming and breaking 
up on the time scale of longer than 20 js at 100 eV (Figure 
S3 in the supplementary [45]). With the moving cluster- 
like ions, "electron bubbles" will transport between dif- 
ferent clusters' interspaces as flowing fluid. These bub- 
bles may crash at one time, but will appear in the next 
time with different configurations in different interspaces. 
The dynamic natures will change rapidly, but the clus- 
ters will persist long enough at the time scale and high 
enough at the percentage to affect the electronic struc- 
tures and related properties such as total energies, op- 
tical properties, and electronic conductivities [43]. This 
understanding violates the traditional assumption that 
the states at so high temperature can be modeled on 
the single atomic scale, and the collective and quantum 
essentials must be treated accurately from first princi- 
ples. Furthermore, the dynamical behaviors of the topo- 
logical structures introduce a challenge for the statistical 
models such as hypernetted-chain (HNC) [16, 17] to in- 
clude more topological networks in its constructed poten- 
tials. From another point of view, the dominant two-ion 
chain structures at the temperature of 100 eV would in- 
duce closer pressures between the statistical methods and 
QMD (QLMD) [25, 26] methods due to the very simple 
structures. 

Let us back to the principal Hugoniot curve in Fig. 4, 
which is determined by interpolating a few density 
points at a fixed temperature according to the Hugoniot- 
Rankine (HR) relation [34, 46]: (/7 - %) = \{P ^ 
^o)(Vb — y\ where /7, P and V are the internal energy, 
the total pressure and the volume of the system, respec- 



tively; f/o, Po, Vb are the respective parameters of the 
initial reference state (7.86 g/cm^ and 20 K). 20 Hugo- 
niot temperature-density points are determined up to the 
pressure of 1 Gbar, as shown in Fig. 4. It is noticeable 
that the spin-polarization is important at T = 0.1 eV. 

The previous data of different experiments are scat- 
tered, which is mainly caused by the uncertainties of the 
temperature and density [5, 28]. Furthermore, different 
models give different results from WDM to HDM, espe- 
cially in the WDM under the pressure of a few Mbar. 
In this regime, our first-principle results are along the 
lower limit envelop of the distributions of the experi- 
mental data, similar to SESAME table. With the in- 
creasing pressures above 10 Mbar, the Hugoniot curve 
derived from QLMD is very close to the experiments 
from rusbank data [27] and within the error bars of 
Batani's experiment [28]. In experiments, external fac- 
tors such as preheating can affect the final results sig- 
nificantly [28], inducing unexpected higher pressures. 
The SESAME tables are far from both experimental 
and our data at relatively high temperature, also shown 
by previous studies [28, 47]. According to the present 
calculations, it can be concluded that below 10 Mbar, 
the experimental results along the lower envelop would 
be regarded error-free, and the limit of pressures in 
experiments above 10 Mbar should locate around our 
points. Similar conclusions can be found in the Hugo- 
niot of Al [34], verifying the accuracy of our calcula- 
tions. The statistical methods such as Thomas- Fermi- 
Dirac (TFD) [26], variational-average-atom-in-quantum- 
plasmas (VAAQP), INFERNO [29] and quotidian equa- 
tion of state (QEOS) [48] models do not give the same ac- 
curate results as QLMD, mainly because physical quan- 
tities such as the energies are sensitive to many-body in- 
teractions and collective quantum electronic distributions 
and their dynamics at high densities and high tempera- 
tures, which are not fully considered in the above statisti- 
cal models. The internal energies from QLMD should be 
lower than those of statistical models according to the 
variational principle. According to HR relation, lower 
internal energies at a definite pressure would result in 
lower densities. Therefore, the Hugoniot curve derived 
from QLMD would be higher than those of the others, 
which is the reason for the results of harder compression 
from first principles [2, 36]. 

Finally, we construct the formula of all calculated EOS, 
i.e., in the range of 0.5 < T < 100 eV, 9 < p < 45 g/cm^. 
This study is important for the applications of our data to 
the experimental benchmark and astrophysical hydrody- 
namics within the wide range of pressures. The optimal 
fit formula in a least-square sense is as follows according 
to the Virial expansions using a typical EOS relation: 

M N 

P{T, P)=Y.CL a^nilogwTDp"', {M = 2,N = 12) 

(1) 
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Where the units of P, T and p are respectively kbar, K 
and g/cm^. The values of the coefficients and the 
validation are shown in the supplementary materials [45] . 

In conclusion, the free electron "blobs" in dense mat- 
ter move and assemble together, forming bigger "bub- 
bles" separated from bond-network clusters at high tem- 
peratures, behaving like quantum flow. The dynamical 
ionic structures are analyzed according to the topological 
structures based on the ECNs idea, giving a new real- 
ization of the stable existence of compact clusters con- 
tributed by the inner-shell electrons forming bonds even 
at T = 100 eV. The ionic structures with bound electrons 
can be compared to the soft "skeleton" in the system, 
which would be responsible for the shearing strength, 
viscosity and response to the shock waves; while the 
free electrons can be compared to "protoplasm" flow- 
ing among the skeletons contributing to kinetic pres- 
sure of electrons, and conductivities. In addition, these 
unique features would raise a challenge to understand 
the dynamical formation of ordered structures from non- 
equilibrium to equilibrium. At last, the new Hugoniot 
data from QLMD simulations give a benchmark on the 
EOS of Fe, which can be regarded as the converged limit 
of the error-free experimental data. 
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